Quantum synchronization of two coupled cavities with second harmonic generation 
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Synchronization is a phenomenon that appears throughout physics, biology, and chemistry. There 
has been much work on how synchronization arises in the classical regime. Motivated by current 
interest in quantum dissipative systems, we investigate whether synchronization can exist in the 
quantum regime. We consider a pair of cavities with second harmonic generation. In the classical 
limit, each cavity has a limit cycle solution, in which the photon number oscillates periodically in 
time. Coupling between the cavities leads to synchronization of the limit cycles. We follow what 
happens to the synchronization as the system becomes more quantum, by decreasing the photon 
number. We find that synchronization between the cavities survives deep in the quantum limit 
when there is much less than one photon in each cavity, because classical correlations are replaced 
by quantum correlations. Our results can be extended to a wide variety of quantum models. 



Introduction. Synchronization refers to the emergence 
of temporal order among a collection of objects, each of 
which oscillate periodically in time |l|-!3!|. Without any 
interaction between them, they would oscillate indepen- 
dently, but the interaction induces them to oscillate in 
unison. Synchronization occurs in many natural settings, 
such as neural networks [4,], audience applause and 
menstration A strikingexample is when thousands 
of fireflies flash in unison [2[. Synchronization also oc- 
curs in man-made settings such as Josephson junctions 
nanomechanical resonators §|, optomechanical ar- 
rays , and trapped ions [ll Jl^ . 

There has been much work on how synchronization 
arises in a group of interacting oscillators. Each mem- 
ber of the group is usually modeled as a nonlinear dy- 
namical system with a limit-cycle solution [ij. A limit 
cycle means that in steady state, the system variables 
oscillate periodically in time. An important feature of 
a limit cycle is that the phase of oscillation is free and 
can take any value, i.e., it is not locked to the phase of 
an external drive. Then due to mutual interaction, the 
limit cycles of different oscillators spontaneously phase- 
lock with each other. In fact, there can be macroscopic 
synchronization in an infinite system, and synchroniza- 
tion can be viewed as a nonequilibrium phase transition 
to temporal order 13l- 15[. Synchronization is inherently 
nonequilibrium, because a limit cycle occurs only in the 
presence of dissipation. 

Models of synchronization are usually based on clas- 
sical equations of motion. Motivated by recent work on 
quantum nonequilibrium systems Il6l4l9l|. we are inter- 
ested in the basic question of whether a quantum sys- 
tem can synchronize. Our approach is to use a quantum 
model that exhibits limit cycles and synchronization in 
the classical limit, and then see what happens as the 
system becomes more quantum. As the quantumness 
increases, one would expect there to be more quantum 
noise [20], which would inhibit synchronization. On the 
other hand, there would also be more entanglement be- 
tween the oscillators, which might enhance synchroniza- 



tion. It is not clear what happens to synchronization 
deep in the quantum limit. 

In this paper, we consider two coupled optical cavi- 
ties, each with second harmonic generation. The classi- 
cal limit here corresponds to many photons in the cav- 
ities. In the classical limit, each cavity has a limit cy- 
cle, in which the photon number oscillates in time (2l|. 
Coupling between the cavities causes the oscillations to 
synchronize in-phase or anti-phase. We study what hap- 
pens as the system becomes more quantum by decreas- 
ing the number of photons in the cavities. As the system 
becomes more quantum, the limit cycles become nois- 
ier. However, classical correlations between the cavities 
are replaced by quantum correlations. We find that syn- 
chronization survives in the extreme quantum limit when 
there is much less than one photon in each cavity. 

We decided to base our work on second harmonic gen- 
eration, because it is the simplest quantum model known 
to have a limit cycle in the classical limit. Our work 
can be extended to the many other quantum models also 
known to have limit cycles, such as Jaynes-Cummings 
cavities ^2, 23], optomechanics ^343], Rydberg atoms 
271 12811 ■ quantum dots [29] , and single-electron transis- 
tors [30|. 

During preparation of this manuscript, we became 
aware of Ref. [3l| . which studies synchronization of op- 
tomechanics in the presence of quantum noise. Our work 
is different, since we are interested in the extreme quan- 
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FIG. 1. Two cavities with nonlinear crystals inside are laser- 
driven and dissipate photons. They are coupled to each other 
due to overlap of their photonic wavefunctions. 
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turn limit, when the classical approximation is completely 
invalid. 

Model. In second harmonic generation, a nonlinear 
crystal within an optical cavity converts light at a funda- 
mental frequency into light at twice the frequency [32 1. 
This setup is widely used for frequency doubling. We ex- 
tend it to the case of two cavities, each with a nonlinear 
crystal inside. The cavities are coupled due to overlap of 
their photonic wavefunctions (Fig. [1} . For simplicity, we 
assume that the cavities are identical. 

We first describe the quantum model and then the clas- 
sical model. Let the two cavities be denoted a and b. 
Each cavity has two modes: ai and 02 are the annihila- 
tion operators for the fundamental and second harmonic 
modes of the first cavity, while hi and &2 are the cor- 
responding operators for the second cavity. An exter- 
nal laser drives the fundamental mode of both cavities, 
and the nonlinear crystals produce light at the second 
harmonic. In the interaction picture and rotating-wave 
approximation, the Hamiltonian is (h = 1) 



H 



iE{a\ 
.X 



ai 



bi) 



al4 + b\%2 - 

+Ai{a\ai + b\bi) + A2(4a2 + 63^2) 
+Vi{a\bi + aibl) + ^2(4^2 + a24). 



(1) 



where E is the laser drive, x is the second-order suscep- 
tibility of the crystal, and Ai and A2 are the frequency 
detunings of the fundamental and second harmonic from 
the laser. For generality and to account for possible 
dispersion, we let Ai and A2 vary independently. The 
term with x means that two photons at the fundamental 
frequency are converted into one photon at the second 
harmonic; the reverse process is also allowed. Vi is the 
coupling between the fundamental modes, and V2 is the 
coupling between the second harmonic modes. 

Since dissipation is a prerequisite for limit cycles, we let 
photons leak out of the cavities. The photon decay rates 
of the fundamental and second harmonic are ki and K2, 
respectively. This open quantum system is described by 
a Lindblad master equation, which says how the density 
matrix p evolves in time: 



P = -i[H,p] + ^ Ki{2aipa\ - a\aip - pa\ai) 

1=1,2 

+ ^ n,{2hpb\~b\hp- pb\b,). 

i=l,2 



(2) 



The master equation_is linear in p and has a unique 
steady state solution 



In the classical approximation to Eq. ([2]) , one assumes 
that each mode is in a coherent state and that there is no 



entanglement between the modes [34|. One replaces the 
operators oi, 02, 61, 62 with complex numbers that denote 
the coherent states, ai, a2, /52- This leads to classical 



equations of motion that are nonlinear: 

E — {ki+ iAi)ai + xo'i(^2 - iVi(3i, 



ai 

Q!2 

$1 



-{k2 + iA2)a2 - |a? - iV2l32, 
-{^2+iA2)P2-\pl-iV2a2. 



(3) 
(4) 
(5) 
(6) 



In the classical model, the average number of photons in 
mode fli is {a\ai) = l^iP and similarly for other modes. 
The classical approximation is an accurate description 
of the quantum model when there are many photons in 
each mode 34] . (This occurs when the laser drive is much 
stronger than the dissipation, since the photon number 
is determined by the balance of driving and dissipation.) 
Intuitively, this is because when a mode is highly popu- 
lated in steady state, it continuously emits photons, so 
an individual photon emission has negligible effect. 

To clarify when the classical model is reliable, it is 
insightful to rewrite Eqs. (IS])-® using scaled variables 
{ai = x^i, Pi = xPi) a-iid scaled parameters {E = Ex)- 
The equations of motion become 



ai 


= E 


- (ki - 


f iAi)<5;i -f 


a*ice2 




(7) 


a2 




-(k2 - 


h 1/^.2)0.2 — 


1-2 


- iV2P2, 


(8) 




= E 


- {^1 - 


f iAi)/3i + 






(9) 


k 




-(k2 - 


h i A2)/32 - 


-Pl~ 
2^1 


iV2a2- 


(10) 



Since Eqs. ©-([TOl) are identical to Eqs. ©-(HI), if we 
vary E and x while keeping Ex constant, the classi- 
cal dynamics are invariant up to a scaling of and Pi . 
This provides a controlled way of following the classical- 
to-quantum transition [3J]. We will solve the quan- 



tum model while decreasing E and increasing x, keep- 
ing Ex fixed, so that the photon numbers decrease, e.g., 
aip ~ jaip/x^. Since the classical dynamics remain 
the same in this procedure, any change in behavior must 
be due to quantum effects. Thus, the classical limit cor- 
responds to large E and small x- The quantum limit 
corresponds to small E and large Xj which makes sense 
since a weak laser drive results in a small photon number. 

To characterize the temporal correlations between the 
various modes, we calculate: 



32(01,02) 



(a|ai4a2) 
(a|ai)(4a2) ' 



{a\a^b%) 

(4«.) 



92{ai,bi) 



(11) 



When 32 > 1, the two modes are positively correlated and 
tend to emit photons simultaneously (bunching). When 
f/2 < 1, they are negatively correlated and tend not to 
emit simultaneously (antibunching). When 172 = 1, there 
are no correlations. 92 is a convenient measure, since 
it works in both quantum and classical models. In the 
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FIG. 2. Classical and quantum trajectories for one cavity, 
showing photon numbers of mode ai (thick, black line) and 
mode 02 (thin, red line) over time, (a) Limit-cycle solution 
of the classical model for E = 8ft:i, x ~ 0.8fi:i. (b)-(d) Quan- 
tum trajectories for the same system, but as it becomes more 
quantum: (b) E = 8ki,x ~ 0.8ki, (c) E = 2ki,x = 3.2ki, 
(d) E — iii,x = 6.4/s:i. All plots use K2 ~ 0.5ki,Ai = 
0.5ki, A2 = Ki. In (d), the modes are antibunched with each 
other due to quantum correlations, despite appearing to be 
positively correlated in the plot. 



quantum model, expectation values are taken with re- 
spect to the steady state p. In the classical model, ex- 
pectation values are averages over time. 

One cavity, classical model. Here we review the results 
of the classical model for one cavity, i.e., Eqs. ^ and ^ 
with Vi = V2 = 0. Suppose one increases E while keep- 
ing the other parameters constant. At a critical value 
Ec, a Hopf bifurcation occurs [2l|. When E < Ec, the 
system has a stable fixed point solution, ai(t) = a\ and 
02(0 = a2, so the photon numbers are constant in steady 
state. When E crosses i?c, the fixed point becomes unsta- 
ble, and a stable limit cycle appears. In this limit cycle 
solution, a\ and ai oscillate periodically in time, and as 
a result, the photon numbers \a\^ and |a2p also oscil- 
late [Fig. [JJa)]. As E increases further, the oscillations 
increase in amplitude and become nonsinusoidal. E^ is a 
complicated function of x, ^i- 

The existence of the limit cycle can be intuitively un- 
derstood as follows. When E is small, ol\ and 02 are 
small, so the nonlinear terms proportional to x in Eqs. (jS]) 
and (HJ have negligible effect. But for sufficiently large 
i?, a\ and a2 are large, and the x terms dominate. The 
effect of the x terms is to exchange energy back and 
forth between the two modes. This exchange is seen in 
Fig. [21 a), where |ai|^ and |a2p are roughly anti-phase 
with each other. 

The temporal correlations between |aip and |a2p are 
captured by (72(01,02)- When E < Ec, there are no cor- 



relations: 32(01,02) = 1. When E > E^, the modes are 
negatively correlated: 32(01,02) < 1. Thus, 32(01,02) is 
an indicator of the existence of the limit cycle. Note that 
the correlations are completely classical. 

Although |aip and ta2p oscillate together, they are 
not said to be "synchronized," since they belong to the 
same limit cycle. Also, it is important to note that the 
limit cycle's phase of oscillation is free and not fixed to 
the phase of the laser drive; this allows two oscillators to 
spontaneously synchronize with each other 

Two cavities, classical model. Now we study the syn- 
chronization of two cavities in the classical limit. Follow- 
ing convention, we assume that the coupling Vi is small, 
so that the limit cycle of each cavity remains well-defined 
and retains its identity in the presence of coupling [l|. 

To identify the steady state solutions of the coupled 
system, it is useful to reduce Eqs. (©-([I]) to the normal 
form of a Hopf bifurcation [s^. This means to expand 
the equations perturbatively in E — Ec and Vi. Again, 
let the fixed point solution of one cavity in the absence of 
coupling be ai and 6:2, which depend on E, x, Hi, Ai. Let 
X and y be deviations each cavity from the fixed point: 



X = 



/ Re(ai - ai) \ 
Im(Q!i — ai) 
Re(Q!2 - 0,2) 

\ Im(a2 - "2) / 



y = 



( Re(/3i - ai) \ 
Im(/3i - ai) 
Re(/32 - 0,2) 

\ Im(/32 - 02) / 



.(12) 



Then x = Za{t)p + Zaityf and y = Zb{t)p + zi,{t)*p*, 
where p is the eigenvector of the linearized system that 
becomes unstable at Ec. Za{t) and Zb{t) are the complex 
amplitudes for each cavity, and they obey 

dz 

= {hi + ih2)Za + {hs + ih4)\Za\'^Za + (/15 + ihG)zb, 

(13) 

(14) 



^ = {hi+ ih2)zb + {h^ + ih4)\zb\'^Zb + (/ig -|- ihe)za 
dt 



where hi depends on the original parameters. 

We are interested in the steady states of Eqs. ([T^ and 
([H]). The two most important steady state solutions are 
Za{t) = -Zb{t) = cie'^i* and Za{t) = Zb{t) = C2e'^^\ 
where ci, C2, f^i, ^^2 are constants. The first solution 
means that |aip oscillates anti-phase with |/3ip, and 
|a2p osciUates anti-phase with |/32p [Fig.[3lja)]. The sec- 
ond solution means that jaip oscillates in-phase with 
|/3ip, and |Q:2p oscillates in-phase with |/32p [Fig. [3l[b)]. 
Which solution is stable depends on the parameters. The 
phase diagram is shown in Fig. SJ^a). 

Equations (fT3)l and (|14p have other steady state solu- 
tions. However, in the limit of small coupling, only the 
two above solutions exist. This is a well-known result 
found by rewriting the equations in terms of amplitudes 
and phases and then adiabatically eliminating the am- 
plitudes to obtain a phase equation [1]. Hence, we limit 
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FIG. 3. Steady state solutions of the classical model for two 
coupled cavities. Photon numbers of mode ai (solid, black 
line) and mode 6i (dashed, red line) are plotted, (a) Anti- 
phase synchrony for Ai — 0.5ki,A2 = ki. (b) In-phase 
synchrony for Ai — — 0.5ki,A2 = — ki. Both plots use: 
E = 8fi:i,x = 0.8ki,K2 = 0.5ki, Vi = 0.1ki,V2 = 0. 
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FIG. 4. Phase diagram for two coupled cavities as a function 
of Ai and A2, showing when ai and fei are bunched or anti- 
bunched, (a) Classical limit, showing in-phase or anti-phase 
synchrony. There is sometimes bistability between the two 
solutions. The parameters are x = 0.8«:i,fi:2 = 0.5fi:i,Vi — 
0.05ki , V2 = 0. For each value of Ai and A2, E is set to Ec + 
O.Sfti. (b) Extreme quantum limit for the same parameters, 
except with S — > and x ^ 00. 



our discussion to the in-phase and anti-phase solutions. 
(The solutions for general coupling have been discussed 
for the cases of ^15 = 12] and he = (36|.) 

Thus, although the coupling is weak, the cavities syn- 
chronize with each other. When the synchronization is 
in-phase, (72(01, ^1), 32(02, ^2) > 1- When the synchro- 
nization is anti-phase, 52(01, ^i), 52(02, ^2) < 1- Again, 
these correlations are classical. 

One cavity, quantum model. Now we consider what 
happens with one cavity in the quantum limit. Figures 
[2I^b)-(d) show simulations of the quantum model, gener- 
ated using the method of quantum trajectories 37, 38| . 
Clearly, as the system becomes more quantum, the limit 
cycle becomes noisier, since the classical approximation 
is less accurate [2^ . (A single photon emission has more 
effect on the wavefunction when there are fewer photons 
overall.) In the extreme quantum limit, when there is 
much less than one photon in each mode, the limit cycle 
is not visually identifiable at all [Fig. [Hd)]. One might 
think that the correlations between ai and a2 have dis- 
appeared. On the contrary, the correlations are strong 
but now quantum in nature. 

To find these quantum correlations in the extreme 



FIG. 5. The correlation (72 ((11,02) for one cavity, plotted as a 
function of E and x, using color scale on right. The black line 
indicates the location of the Hopf bifurcation Ec as a function 
of X- Parameters are K2 = 0.5ki, Ai = 0.5ki, A2 = ki. This 
plot was generated using the quantum trajectory method, so 
there is some noise due to statistical uncertainty. 



quantum limit, we analytically solve for the steady state 
p of Eq. ([2]) in a perturbation series of E. By doing this 
to sixth order in E, we find 52(01, 02) to first order in E: 

52(01,02) = l.,4^,.r A./A.^A.>^......^...M +0{E^). 



-, I ix^+X^[-Ai(Ai+A2)-|-Bi(>ii-l-K2)] 

(Af + Kf)[(Al+A2)2+(Kl+K2)2] 



(15) 



Recall that the quantum limit corresponds to i? — ?> 
and X — )■ 00. Thus, in the extreme quantum limit, 
52(01,02) = 0, i.e., there is strong antibunching between 
ai and 02- This result is independent of Ki, A^. 

Figure [5] shows the classical-to-quantum transition. In 
the classical limit, 52(01,02) = 1 when E < Ec and 
52(01,02) < 1 when E > Ec, where Ec is the location 
of the Hopf bifurcation. As the system becomes more 
quantum, the transition at Ec smoothes out, and the 
region of antibunching expands into E < Ec- Also, an- 
tibunching in the classical limit gradually develops into 
strong antibunching in the quantum limit. In this sense, 
the two modes still exchange energy back and forth, but 
the correlations are now quantum instead of classical. In 
the extreme quantum limit, the classical approximation 
is completely wrong, demonstrating the presence of quan- 
tum correlations. In fact, the presence of entanglement 
is clearly seen in Fig. EJd): when one mode decays, the 
other one instantaneously changes. 

Two cavities, quantum model. Finally, we consider the 
synchronization of two cavities in the quantum limit. Re- 
call that in the classical limit, the cavities were in-phase 
[52(01,61) > 1] or anti-phase [52(01,61) < 1] with each 
other, depending on the parameters. We want to see 
what happens to these correlations in the extreme quan- 
tum limit. In classical systems, noise inhibits synchro- 
nization because it causes the phases of the oscillators 
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to drift with respect to each other Hence, one would 
expect quantum noise to wash away any trace of synchro- 
nization between the oscillators. 

To find 172(01, bi) in the quantum limit, we again solve 
Eq. ([2]) in a perturbation series of E but now include Vi 
and V2 ■ By solving for the steady state p to fourth order 
in E, we obtain an analytical expression for (72(01, ^1) to 
first order in E. In the limit of x — >■ 00, 

92iaM = ^^^^^^^ + 0{E% (16) 

which is independent of A2, K2, V2. For small Vi, 

2A1F1 



52(01,61) ~ 1 



(17) 



Whether oi and bi are bunched or antibunched depends 
on the sign of AiVl. The phase diagram is shown in 
Fig. mjb). (Numerically, we find that the phase diagram 
for a2 and 62 is very similar.) The fact that the phase 
diagram is different from the classical result in Fig. |4lja) 
demonstrates the importance of quantum correlations in 
the quantum limit. 

Equation ()17p shows that even in the extreme quantum 
limit, there are significant correlations between ai and 
bi. So despite the presence of much quantum noise, syn- 
chronization survives because classical correlations are 
replaced by quantum correlations. Note that while one 
might expect some sort of correlation for finite E and x, 
it is not obvious that the correlations would still exist in 
the limit of i? — > and x 00. 

Conclusion. We have shown that synchronization sur- 
vives deep in the quantum limit. Our result raises the 
question of whether a macroscopic number of oscillators 
can also synchronize in the quantum limit. Classical syn- 
chronizing systems are known to exhibit phase transi- 
tions and scaling behavior similar to equilibrium systems 
[39I [io! ]. One should investigate how quantum fluctua- 
tions affect these critical properties. 

It would be interesting to study quantum synchroniza- 
tion using the other quantum models known to exhibit 
limit cycles, as discussed in the introduction. A promis- 
ing approach is to use coupled cavities with Jaynes- 
Cummings nonlinearity 18|, |41| , since one can fabricate 
a large of array of such cavities in the microwave regime 
[4^ . Perhaps synchronization is as widespread in the 
quantum regime as it is in the classical regime. 
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